Supernova-driven outflows and chemical evolution 
of dwarf spheroidal galaxies 

Yong-Zhong Qian * and Gerald J. Wasserburg ^ 

*School of Physics and Astronomy, University of Minnesota, Minneapolis, MN 55455, and ^The Lunatic Asylum, Division of Geological and Planetary Sciences, California 
Institute of Technology, Pasadena, CA 91125 

Submitted to Proceedings of the National Academy of Sciences of the United States of America 



We present a general phenomenological model for the metallicity 
distribution (IVID) in terms of [Fe/H] for dwarf spheroidal galaxies 
(dSphs). These galaxies appear to have stopped accreting gas from 
the intergalactic medium and are fossilized systems with their stars 
undergoing slow internal evolution. For a wide variety of infall his- 
tories of unprocessed baryonic matter to feed star formation, most 
of the observed MDs can be well described by our model. The key 
requirement is that the fraction of the gas mass lost by supernova- 
driven outflows is close to unity. This model also predicts a rela- 
tionship between the total stellar mass and the mean metallicity for 
dSphs in accord with properties of their dark matter halos. The 
model further predicts as a natural consequence that the abundance 
ratios [E/Fe] for elements such as O, Mg, and Si decrease for stellar 
populations at the higher end of the [Fe/H] range in a dSph. We 
show that for infall rates far below the net rate of gas loss to star 
formation and outflows, the MD in our model is very sharply peaked 
at one [Fe/H] value, similar to what is observed in most globular 
clusters. This suggests that globular clusters may be end members 
of the same family as dSphs. 
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In this paper we show that supernova-driven gas outflows 
play a prominent role in the chemical evolution of dwarf 
spheroidal galaxies (dSphs) . In the framework of hierarchical 
structure formation based on the cold dark matter cosmology, 
dwarf galaxies are the building blocks of large galaxies such 
as the Milky Way. In support of this picture, some recent ob- 
servations showed that elemental abundances in dSphs of the 
Local Group match those in the Milky Way halo at low metal- 
licities (e.g., [1, 2, 3, 4, 5]; see [6] for a review of earlier works). 
It is expected that detailed studies of chemical evolution of 
dwarf galaxies can shed important light on the formation and 
evolution of the Milky Way in particular and large galaxies 
in general. Here we present an analysis of the evolution of 
[Fe/H] = log(Fe/H) - log(Fe/H)0 focussing on dSphs. The 
approach is a phenomenological one that takes into account 
infall of gas into the dark matter halos associated with these 
galaxies, star formation (SF) within the accumulated gas, and 
outflows driven by supernova (SN) explosions. The sources for 
production of Fe are core-collapse SNe (CCSNe) from progen- 
itors of 8-100 A'Iq and Type la SNe (SNe la) associated with 
stars of lower masses in binaries. Observations require that 
some SNe la must form early along with CCSNe without a 
significant delay. It will be shown that there is a direct and 
simple connection between the metallicity distribution (MD) 
for a given dSph and two parameters Apo/A and a. The ratio 
refers to the net rate Ape of Fe production and the net rate A 
of gas loss to SF and SN-driven outflows, and a indicates the 
promptness for reaching peak infall rates. This model explic- 
itly predicts the ratio of the stellar mass in the dSph to the 
total mass of the host dark matter halo. 

Phenomenological models for chemical evolution have a 
long history (e.g., [7]) and were applied to dSphs previously 
(e.g., [8, 9, 10]). Dynamic models for dSphs including dark 
matter were also studied (e.g., [11, 12]). The first effort was 
made in [13] to reconcile models of hierarchical structure for- 



mation involving dark matter halos with the then-available 
luminosity-radius-metallicity relationships for dwarf galaxies. 
There it was shown that SN-driven outflows could explain the 
observed trends. Recent observations [14] give results for eight 
dSphs with rather detailed structure of their MDs and pro- 
vide a basis for exploring models of their chemical evolution 
(e.g., [15]). Among the key issues that we try to address are 
MDs exhibited by stellar populations of dSphs. Our general 
approach follows that of Lynden-Bell described in his incisive 
and excellent article on "theories" of the chemical evolution 
of galaxies [16] . It will be shown that the metallicity at which 
the MD peaks is directly related to the efficiency of SN-driven 
outflows and that the MD of a dSph and the relationship 
between the stellar mass and the mean metallicity for these 
galaxies are direct consequences of the model. These results 
are in strong support of those of [17], where earlier and less 
precise data on metallicities of dwarf galaxies were used to 
address this problem. 

In our approach, we consider evolution of Fe in a homo- 
geneous system of condensed gas governed by 
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where Mg{i) is the mass of gas in the system at time t, 
{dMg / dt)\n is the infall rate of pristine gas, ■i/>(f) is the star for- 
mation rate (SFR), -F'out(t) is the rate of gas outflow, Mpo(t) 
is the mass of Fe in the gas, and PFc{t) is the net rate of Fe 
production by all sources in the system. We assume that the 
SFR is proportional to the mass of gas in the system with an 
astration rate constant A,, 
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Given {dMg/dt)in, Fout{t), and Ppcit), Eqs. 1 and 2 can be 
solved with the initial conditions Mg{0) = and Afpc(O) — 0. 

The MD of a system measures the numbers of stars formed 
in different metallicity intervals that survive until the present 
time. We use [Fe/H] to measure metallicity. As the mass frac- 
tion of H changes very little over the history of the universe, 
we take [Fe/H] — logZpc, where 
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Here Xjj^ is the mass fraction of Fc in the sun. We assume 
that the initial mass function of SF does not change with time 
and is of the Salpeter form. Then the number of stars formed 
per unit mass interval per unit time is related to the SFR as 
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where m is the stellar mass in units of Mq with mi and niu be- 
ing the lower and upper limits, respectively. We take mi = 0.1 
and mu = 100. Assuming that Zpeit) increases monotonically 
with time, we obtain the MD 
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where mmax(t) is the maximum mass of those stars formed at 

time t that survive until the present time. There is little SF 
in dSphs at the present time (cf. [18]). We assume that SF 
ended at time t ; in a system. Then the total number of stars 
in the system at the present time is 
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The integral involving mmax(i) in Eqs. 6 and 7 increases only 
by 6% when mniax(t) increases from 0.8 to 100. As stars with 
m — 0.8 have a lifetime approximately equal to the age of the 
universe, we take mmax(t) = 0.8 in both these equations to 
obtain the normalized MD 
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The Model 

The key input for our model is the infall rate {dMg/dt)in, the 
outflow rate Fout(t), and the net Fe production rate Ppeit). 
The latter two rates are closely related to the occurrences of 
CCSNe and SNe la in a system. CCSNc arc associated with 
massive stars (8 < m < 100) that evolve rapidly. In con- 
trast, SNe la require consideration of the evolution of binaries 
involving lower-mass stars with longer lifetimes. In all our 
previous studies (e.g., [19]), we considered that the evolution 
timescale for SNe la was ~ 1 Gyr using the lifetime of stars 
with m ~ 2. This meant that SNe la would not contribute 
to the Fe inventory during early epochs and was in accord 
with the general approach used by other workers. The conse- 
quence of this assumption is that the MD for a system must 
have two peaks due to the assumed late onset of SNe la (e.g., 
[19]). However, of the eight dSphs studied in [14], only a single 
peak is observed in the MD for Fornax, Leo I, Leo II, Sextans, 
Draco, and Canes Venatici I, and there is only some indicar 
tion for two peaks in the MD for Ursa Minor and perhaps 
Sculptor. Prom this we conclude that there must be a prompt 
component of SNe la that start in a system on much shorter 
timcscales than ~ 1 Gyr. Such a component is supported by 
both supernova surveys (e.g., [20]) and models that consid- 
ered detailed evolution of various binary configurations (e.g.. 



[21, 22]). The occurrences of "prompt" and "delayed" SNe la 
were investigated in [23], where it was argued that both pop- 
ulations were present at high redshift. The nature of these 
two classes of SNe la remains unclear. For simplicity, we will 
lump the Fe production by these sources together with that 
by CCSNe and ignore the time delay between the birth and 
death of all SN progenitors, except when the amount of gas in 
the system is so low that SNe la would dominate. A detailed 
treatment of SNe la that takes into account their evolution 
timescales will be discussed in a subsequent paper. 

Under our assumption, the total rate of CCSNe and SNe 
la, and hence JVe(i), are proportional to the SFR. We further 
assume that the rate of outflows driven by these SNe is also 
proportional to the SFR. Specifically, we take 
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where is a dimensionless constant that measures the effi- 
ciency of the SN-driven outflows and Ape is a rate constant 
that is proportional to A* and the effective Fe yield of SNe. 
We take the infall rate to be 
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where Ain is a rate constant, Mq is the total mass of gas infall 
over < t < GO, and T{a + 1) with a > — 1 is the Gamma 
function of argument a -|- 1. The above modified exponential 

form was specifically chosen to explore the role of the time 
dependence of the infall rate in chemical evolution (cf. [7]). 
The infall rate peaks at t = for —1 < a < 0, and the peak 
time increases to a/Ain for a > 0. The form with a > allows 
a slow start of significant gas accumulation in the system. 
With the above assumptions, Eqs. 1 and 2 become 
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where A = (1 + 'q)\^. Note that for Ain <C A, the solutions 
to the above equations approach a secular state for which 
Mgit) ~ {dMg/dt)in/X and Zpe ~ Ape/A. This is analogous 
to the quasi-steady state for the metallicity of the interstellar 
medium first proposed in [24]. For simplicity, we assume that 
Ain and A are so large that tf — oo can be used effectively 
in Eq. 8 for the MD. As Mg(0) = Mg(oo) = 0, integrating 
Eq. 12 over t gives A Mg{t)dt = Mq. So the total gas mass 
used in SF is A* Mg {t)dt = Mo /(I + r;) and the remainder 
of the gas infall is blown out as outflows. We assume that all 
outflows are lost from the system into the broader intergalae- 
tic medium (IGM), thus enriching the latter in metals. This 
approach gives a natural cutoff to the chemical evolution of 
the system when the total mass of gas lost to SF and out- 
flows is equal to the total mass of gas infall (partial recycling 
of outflows would increase the degree of chemical enrichment 
but is ignored here for simplicity). At the present time, the 
total mass of stars in the system can be estimated as 
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where we have used Mo = {Q,b/Q.m)Mh = O.ITM^. Here Mh 
is the total mass inside the dark matter halo hosting the sys- 
tem, and Q,b and 0„i are the fractional contributions to the 
critical density of the universe from baryonic and all matter, 
respectively. 
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Exploration of Results from the Model 

The parameters governing the solutions to Eqs. 12 and 13 are 
a, Ain, A, and Ape. We illustrate the dependences of the MD 
on these parameters in the following subsections. 

Dependence of the MD on Ain/A with OC — O.For a = 0, 

the solutions to Eqs. 12 and 13 are 
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For Ain /A = 1, the above two equations reduce to 
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In the limit Ain = oo, which is equivalent to setting Mg(0) 
Mo and {dMg / dt)\T, = for t > 0, the results are 
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In general, Eqs. 17 and 18 give the MD as a function 
of [Fe/H] in parametric form. This MD only depends on 
Ain/A and Apo/A, but not on the absolute values of these 
rates (this is true so long as At/ ^ 1 and Aint/ 3> 1, see 
Eqs. 8 and 15). For a fixed Ain/A, changing Apc/A only trans- 
lates the MD along the [Fe/H]-a:x;is. This can be most easily 
seen in the special cases of Ain/A = 1 and Ain = oo, where 
N^JdN/dlFe/B] is a simple function of At and [Fe/H] differs 
from log(At) only by a shift of log(Ape/2A) and log(Ape/A), 
respectively. The shape of the MD is determined by Ain/A. 
This can be seen from Fig. lA, which uses Apc/A — 0.1 
and shows that as Ain/A increases from 1/2 to oo, the po- 
sition of the peak of the MD changes slightly [but staying 
close to [Fe/H] — log(Ape/A) = —1, the exact peak posi- 
tion for Ain/A = 1 and Ain = oo] and the shape of the MD 
becomes broader. Note the sharp cutoff of the MD to the 
right of the peak for Ain/A = 1/2 with no stars formed above 
[Fe/H] = log[Ape/(A - Ain)] = -0.7. 

The case of Ain/A < 1/2 requires separate discussion. For 
illustration, we again take Apo/A — 0.1 and show the MD 
for Ain/A = 0.1 in Fig. IB. This MD has an extremely sharp 
peak with 90% of the stars having -1.154 < [Fe/H] < -0.954. 
This is simply a case close to secular equilibrium (see discus- 
sion below Eqs. 12 and 13; cf. [24]). Piling up of stars in an 
extremely narrow metallicity range is typically not observed 
for dwarf galaxies but strongly resembles what is observed for 
most globular clusters. In general, for Ain/A significantly be- 
low 1/2, stars are concentrated immediately below [Fc/H] = 
log[Ape/(A — Ain)]. For comparison, the solid curve in Fig. lA 
shows that the MD for Ain/A = 1/2 first rises to a peak and 
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then sharply drops to zero as [Fe/H] — ;> log[Ape/(A — Ain)]. 
This general behavior also applies to Ain/A > 1/2 but with a 
more extended tail at high metallicities for a larger Ain/A (see 
Fig. lA). 
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Fig. 1. Example MDs for Ape/A = 0.1 and various values of Ain/A with 
a = 0. Note that in general the shape of the MD is determined by Ain/A and chang- 
ing Ape/A only translates the MD along the [Fe/H]-axis. (A) The solid, dashed, dot- 
ted, and dot-dashed curves are the MDs for Ain/ A = 1/2, 1, 2, and oo, respectively. 
The vertical dashed line indicates the peak at [Fe/H] = log(AFe/A) = —1 for the 
dashed and dot-dashed curves. (B) The solid curve is the MD for Ain/A = 0.1 and 
the dashed line indicates the limiting value of [Fe/H] = log[AFe/(A — Ain)] = 
—0.954. Concentration of stars in an extremely narrow range of [Fe/H] as shown 
by the solid curve is typically not observed for dSphs but strongly resembles what is 
observed for most globular clusters. 



Dependence of the MD on a with Ain — A. Based on the 
above discussion and the example MDs shown in Fig. 1, it 
seems reasonable to choose Ain = A and explore possible MDs 
for different values of a and Ape/A. For Ain = A and a > —1, 
the solutions to Eqs. 12 and 13 are 
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which give 
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The above MD again has a peak at [Fe/H] = log(AFe/A) for 
all a > —1. Changing Afc/A only translates the MD, shifting 
the peak in particular, along the [Fe/H]-axis. The shape of 
the MD is determined by a. For a specific Apo/A, the MD 
becomes narrower while peaking at the same [Fe/H] when a 
increases above —1. A positive a reflects slower infall at the 
start of the system, which suppresses SF at the early stages 
and reduces the extent of the low-metallicity tail of the MD. 
For a < 0, the initial infall rate is enhanced and consequently, 
more stars are formed at lower metallicities. 

Note that if we take the limit a = —1, the MD given by 
Eqs. 25 and 26 coincides with that for Ain = oo and a — (see 
Eqs. 21 and 22). Thus, for Apc/A = 0.1, as a increases from 
— 1 to 0, the MD for Ain ~ A changes from the dot-dashed to 
the dashed curve shown in Fig. lA while peaking at the same 
[Fe/H] = — 1. More example MDs for Ain ~ A and a > —1 are 
shown in Fig. 2 where we compare them with observations of 
dSphs. 

Comparison with Observations 

We now explore the implications of our model for observations 
of dSphs. We first discuss the MDs using the high-quality data 
set of [14] and then study the relationship between the mean 
metallicity ([Fe/H]) and the stellar mass M, of dSphs. 

MDs for dSphs. Important medium-resolution data on the 
MDs for dSphs were provided in [14] and are summarized in 
Fig. 2. Careful inspection of these data shows that there is 
only some indication for two peaks in the MD for Ursa Minor 
and perhaps Sculptor. An MD with two peaks would be typi- 
cal if the turn-on of SNe la were sudden and with a significant 
delay relative to CCSNe. It was the lack of such MDs from 
observations that led us to pursue a model where the net Fe 
production rate is without discontinuities. 

We focus on models with Ain = A and different values of 
a. In this case, a model MD is specified by Ape/A and a, 
which determine its peak position [Fe/H] = log(AFe/A) and 
its shape, respectively. Positive values of a corresponding to 
slower initial infall give rise to narrower MDs while negative 
values corresponding to more rapid initial infall result in more 
extended MDs. As the MD is normalized, its peak height can 
be used to estimate a effectively. Using the position and the 
height of the peak for the observed MD as guides, we fit a 
model MD for each of the eight dSphs reported in [14]. The 
observed and fitted MDs are shown as histograms and curves, 
respectively, in Fig. 2. The adopted values of Ape/A and a 
are indicated for each galaxy. These values were not obtained 
from the best fits, but were simply picked to illustrate the 
overall adequacy of our model. Very good fits are obtained 
for Fornax and Leo I, which arc the most massive of the eight 
dSphs. The fits for the least massive four, Draco, Sextans, 
Ursa Minor, and Canes Venatici I, are rather good, although 
the model MDs appear to underestimate their stellar popula- 
tions at the highest metallicities. The only exceptions are Leo 
II and Sculptor with intermediate M* . Using the same Ape/A, 
we obtained a better fit to the data for Leo II with Ain/A = 0.8 
and a = (dashed curve in Fig. 2C) than with Ain/A = 1 and 
a = 1 (solid curve) . However, neither type of MD can provide 
a good fit to the data for Sculptor, which may represent an 
MD with two peaks. There is perhaps a more clear indicar- 



tion for such an MD for Ursa Minor. We note that if infall 
took a more complicated form than Eq. 11, then the shape 
of the MD would change accordingly. It is possible that the 
broadened peak in the MD for Sculptor might be explained 
by an increase in the infall rate after the assumed smooth rate 
in Eq. 11 peaked. However, in no case can a significant part 
of the enriched outflows be returned to the infalling matter 
as this would produce a population of stars with high [Fe/H] 
values, which is not observed. 

In a previous study [14] with alternative modelling, de- 
tailed multi-parameter fitting was used. As can be seen from 
Fig. 2, a good to excellent fit of the model to the data can be 
obtained in our approach. The model thus appears to give a 
good description of the MDs of dSphs and is easily understood 
in terms of physical processes. 

Relationship between ([Fe/H]) and M* for dSphs. The ob- 
served MDs of dSphs require that Ape/A ~ 0.01-0.1 (see 
Fig. 2). As Ape/A ~ (1 -|- 77)"^ (see derivation of Ape in 
SI Text), this requires that rj ~ 10-100. Consequently, only a 
fraction {l+ri)~^ 1-10% of the gas falling into the dark mat- 
ter halo hosting a dSph is used in SF and the rest is blown out 
of the halo by SN-driven outfiows. The MD of a dSpli peaks 
at essentially its mean metallicity ([Fe/H]) ~ log(Apc/A) ^ 
— log 77. The lower ([Fe/H]) a dSph has, the lower fraction 
of the infalling gas is stored in its stars. There is thus a di- 
rect relationship between ([Fe/H]) and the stellar mass M, 
for dSphs. Using the data on ([Fe/H]) in [14] and those on 
M, in [25, 26], we show this relationship in Fig. 2. The dot- 
dashed line in Fig. 2 has a slope of 2.5 and passes through the 
point defined by the average values of the data. It is nearly 
the same as the least-square fit (solid line) and is in excellent 
agreement with the result in [25] (sec SI Text). 

As discussed in [17], the slope of 2.5 for the relationship 
between log(Mt/M0) and ([Fe/H]) has a simple physical ex- 
planation. In general, this slope follows from the relationship 
between the radius rh and the total mass Mh of the dark 
matter halo hosting a dSph in addition to the dependences of 
M, and ([Fe/H]) on rj discussed above. The outflow efficiency 
rj is inversely proportional to the depth of the gravitational 
potential well of the dark matter halo. 
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where we have assumed rn oc in the second step. This 
gives 



([Fe/H]) ~ - log 7? ~ (1 - /3) log Mh + const. 
In addition, Eq. 14 gives 
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Thus, a slope of 2.5 for the relationship between log(M,/M0) 
and ([Fe/H]) requires (2-^)/(l-/3) = 2.5, or ^3 = 1/3, which 
is in agreement with the framework of hierarchical structure 
formation (e.g., Eq. 24 in [27]). 



Conclusions 

We have shown that the MDs of dSphs can be reasonably well 
described quantitatively by a phenomenological model based 
on general considerations of gas infall into a dark matter halo, 
astration in the condensed gas, and SN-driven gas outflows. 
The peak position of the MD is governed by the ratio of the 
rate constants for net Fe production and net gas loss to astra- 
tion and outflows, Ape/A = Ape/[(1 -|- 7?)A«]. The observations 
require high efficiency rj ~ 10-360 (see Table SI) for SN-driven 
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Fig. 2. Comparison of model MDs with observations of dSphs. The data are taken from [14] and shown as histograms with error bars. The model MDs assume the indicated 
parameters and are shown as curves. The dashed curve in (C) provides a better fit to the data than the solid curve. Values of M* are taken from [25] for (A)-(G) and from 
[26] for (H). Note that if the baryonic matter is not always blown out of the dark matter halo but returns to the gas mass in a dSph after some time, then the curves will have 
a leading-edge tail going to higher [Fe/H]. 



outflows and hence, massive gas loss from dark matter halos 
associated with dSphs. The model also directly relates the 
stellar mass A/, remaining in a dSph to its mean mctallicity 
([Fe/H]) through the efficiency of outflows governed by the 
mass Mh and the radius rh of the dark matter halo. The 
observed relationship between log(M,/MQ) and ([Fe/H]) has 
a slope of 2.5 over the wide range of 4.8 X 10^ < M./Mq < 
4.6 X 10* for dwarf galaxies (sec SI Text). This indicates 

fh oc M^^"* in agreement with the framework of hierarchical 
structure formation. Our results confirm the previous studies 
of [17] and [25], and are in support of the early work of [13] 
on the general model of dwarf galaxy formation. 

Our model also demonstrates that for slow infall rates 
Ain <^ A, the resulting MD must be extremely sharply peaked 
(essentially concentrated at a single value of [Fe/H]). Such 
MDs apply to all globular clusters with the exception of the 
most massive ones (e.g., [28]). This suggests that globular 
clusters are the result of the general process of astration gov- 
erned by very slow infall rates compared to net gas loss rates. 
Those slow infall rates further suggest that globular clusters 
might be weak feeders during the inliomogcncous evolution 
inside a large system and that they are not responsible for sig- 



nificant depletion of the baryonic supply of that system. As 
such, the dark matter halo would be associated with the entire 
system, but not bo specifically tied to the globular clusters. 
This is in contrast to dSphs, which result from processing all 
the baryonic supply in their dark matter halos and are natu- 
rally associated with these halos. Nonetheless, from the point 
of view of the analysis presented here, we consider globular 
clusters to be part of a family of early-formed dwarf galaxies 
but with low effective infall rates compared to outflow rates. 
With regard to the diverse morphological types of dwarf galax- 
ies, we consider dSphs to represent isolated evolution without 
dynamic effects from mergers or tidal interactions with nearby 
systems. Other morphological types may result from such ef- 
fects. Insofar as all dwarf galaxies are related by the same 
general process as presented here, the observed ongoing as- 
tration in some dwarf irregular galaxies (e.g., [18]) poses a 
problem. These systems must be experiencing secondary pro- 
cesses of gas accretion due to local infall or mergers. Their 
MDs should have a second peak due to the late infall. 

The analysis presented here assumes that there is not a 
discontinuous onset of SNe la contributing Fe. Otherwise, the 
general outcome would be MDs with two peaks, which are at 
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most only rarely observed. It is not appropriate to extend 
this analysis to other metals than Fc without a more real- 
istic treatment of the relative contributions of CCSNc and 
SNe la that each produce very different yields of the other 
metals. Such a treatment will be carried out in a subse- 
quent paper that takes into account the detailed distribution 
of the delay between the birth and death of the progenitors 
for SNe la. We note here that as the rate of CCSNe decreases 
with decreasing gas mass, the contributions to Fe from SNe la 
will become larger or dominant. This is because previously- 
formed stars of low to intermediate masses in binaries will 
continue to evolve and produce SNe la even after the SFR 
decreases. Crudely speaking, the net Fo production rate ap- 
propriate for later times must be changed from Eq. 10 to 
PFe(t) ~ X® [Agf M<,(t) + \l%Mg{t - A)], where X^^ and A^.^ 
are the rate constants for CCSNe and SNe la, respectively, 
and A is a typical delay between the birth and death of the 
progenitors for SNe la. As SNe la do not produce e.g., O, Mg, 
and Si, this will result in lower values of [E/Fe] for these ele- 
ments at higher values of [Fe/H], which is in agreement with 
the general trend observed in dSphs (e.g., [15]). 

In conclusion, dSphs have evolved as a result of very mas- 
sive gas loss and this gas has gone into the medium outside the 
dark matter halos associated with these galaxies. This means 
that these dark matter halos must have masses ~ 10^-10^ 
times greater than the present total mass in stars after the gas 
loss and the cosmic ratio of baryonic to all matter are taken 
into account (see Table SI). This is in qualitative accord with 
the conclusions from previous studies to infer the halo masses 
from observations of dSphs (e.g., [29]). As a result of the mas- 
sive gas loss, extensive enrichment of metals has occurred in 
the general hierarchical structures outside individual dSphs. 
This drastically alters the subsequent chemical evolution of 
the emerging larger galaxies as most baryonic matter must 
have passed through processing in dSphs. As more massive 
dark matter halos axe formed during hierarchical growth, the 
efficiency of outflows will decrease. Nevertheless, the general 



IGM must have been enriched by the net outflows from all 
dSphs. If the average IGM has [Fc/H] ~ —3 and the outflows 
from dSphs have [Fe/H] ~ —1.5 on average, this would imply 
that ~ 3% of all baryonic matter was processed in dSphs. 
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Fig. 3. The relationship between log(M,/Af0) and ([Fc/H]) for dSphs. Val- 
ues of for filled and open diamonds are taken from [25] and [26], respectively. All 
values of ([Fe/H]) used here are taken from [14]. The solid line is the least-square 
fit to the data and the dashed lines indicate the la error in Iog(Mt/M0). The 
dot-dashed line has a slope of 2.5 and is nearly the same as the solid line. 
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SI Text 

RelaticHiship bQtween ([Fe/H]} and A/* fordSphs. Tbedaia 
on {[Fe/H]) and M, for ihe eight dSphs shewn in Fig. 2 plus 
another six from literature siq given in Table SI and shown in 
Fie. 3' AH values of ([Fe/H]) are taken from [1]. The vaiiies 
of AJ, for FornaK, Leo 1, Leo IL Sculptor, Draco, Sex tans, and 
Llrsa Minor ^fifled diamonds in Fig. 3) are t.aken from [2] wht!e 
those for Canes Venatici I, Hercules, Ursa Major [, Leo TV. 
Canes Venatici II, Uisa Major [I, and Coma Berenices (open 
diainoEid&} are taken ffom [3]. (An important input to infer 
the Mt of a dwarf galaxy is the present mass distribution of 
its ste4]ar pcpulaEiona. TVo different mas distribations were 
used in [3], which ga^-e values of M, differing by a factor of 
fB 2. For consistency, we have adopted the results derived 
in [3] from the mass distribution that is similaT to the one 
used in [2].) These values of A/, arc inferred from obscrva- 
tiona and havis some considerable uncertainties, For CKample, 
Draco is common to the data sets of [2] and [3] and is inferred 
to have Aft = 0. 1 x lO'^ Aft and 3.'2 x lO'^ Afg,, rspectively. In 
view of this, we treat laE(A/, /Af a) as a function of the better 
determined {[Fe/H]} and obtain a least-square linear fit 

bE(.lf»/Af£;) = lO.OS ±0.53+ (2.41 ±0.29)([Fe/H]). Jl] 

The sohd Jine in Fig, 3 corresponds to the central vdues of 
the above fit. We estimate that the Icr error in the inferred 
values of log(Af, /Mgj) is 0.47 dex. This error is indicated by 
the mo dashed lines in Pig. 3. 

Tlie relationship in Eq. SI is essentially the same bs that 
foimd in [2], which used a set of more massive dwarf galax- 
ies (4.0 x IQ^^ < Af,/Af(i, < 4.6 X 10^, inclndiELg both ir- 
regular and spheroidal galaxies) and earlier and less precise 
data on {[Fe/H]). That work also used a different parameter 
logZ = {[Fe/H]) -I- log Zii, where Zg, = 0,019 is the total mass 
fraction of metals in the sun [4], In the repr^sntation used 
here, [2] obtained 



loE(Af. /Afa } = 10,37 + 2.S{[Fe/H]). 



[2] 



Our result in Eq. SI for a different mass range (4.8 x 10* < 
Af, /A/s < 1.9 X lO'') using only dSphs is in excellent agree- 
menl with Eq. S2. The dol^dashed line in F^g. 3 has a slope 
of 3.5 and passes through the point defined by the average 
values of the dat,a. It corresponds to 



log(Af,/Af^i) = 10,2e -I- 2.S{[Fe/H]), 



[3] 



and is nearly the same as the solid hne. The mass ranges cov- 
ered by Eqs, S2 and S3 overlap and both include the seven 
dSphs shown as filled diamonds in Fig. 3. Tbgether these two 
results show that log(Af,/A/,ii ) increases with ([Ffe/H]) with 
a slope of 2.0 over the wide range of 4,8 x 10* < Af^/Atfo < 
46 X 10^ for dwarf gataxies. 

Parameters of the Model. We can estimate the parameters of 
our model based on piopenios of the dark matter halo hostiDg 



3 dSph and characteristics of CCSNe and 5Ns la. The rata of 
CCSNe is 



,■100 JlKJ 



Af,(t) 



Afa 



[4] 



For simplicity, we asume that the rate of SNe la is fllia{£) = 
kR(zc{t), where A; is a constant. The average Fe yield of each 
CGSN is 



.-100 v,CC f„^„-1.3B J„ 



[EJ 



where Yg^{7ii) is the Fe yield of a CCSN from a progienitor 
of mass m and we have taken Ypip[TTi) = 3 x lO"' Af^i for 
e < m < 11 IS] and ? X 10"^ Aff. for 11 < m < 100 (eg., 
in the second step. We take the average Fe yield of each SN 
la to be (ypa) = D.J'Afa (a^g., ]T].). The net Fe production 
rate is 

PF.{t) = Ap^^JWs(i ) = (r^^^ )Rcc{ t) + {Yi^} Rl. (£) . [8] 

Using the above results and — 1.3 x 1D~^ [4], we obtain 

Af,j =0.301(1 -|-lS.3fc)A,. [T] 

As both CCSJ^Je Md S!*Je la produce Fe, the assumption that 
Ria.{t)/Rcc{t] is a constant k does not critically affect the 
model presented here for Fe in consideration of the approKi- 
mations made. I-ltHPever, other elements such as O. Mg, and 
Si are produced by CCSNe but not by SNe la. For these ele- 
ments the detailed distribution of the delay between the birth 
and death of the pTogenitors for SNe la must be taken into 
account. This will be dealt with in a subsequent paper that 
treats a number of elements including Mg and Fe with detailed 
considETations of SNe la. 

The rate of SN-driven outflows iaF'iHii(t) — i^AtAfjIt) and 
can be estimated by 



GAf, 



= (£S^)/Ioc(f ) + {£iL)fli.(t), [S] 



Reserved for Pubiicatioii Footnotes 
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where G ia the ^svitaxkmal constant, and (E^n) i^n) 
are Lhe aver^ kiiiQtic energy imparted to thg suTrounding 
each CCSi^ and SN la, respectively. Specifically, 

{£^3.).= TTT^ =6.do>;10 erg, [fi] 

where ££^(j7i) ia the kinetic energy imparted to the surround- 
ing gas by a CCSN from a progenitor of mass m, and we have 
taken E^^{m) = lO*" erg for 8 < m < 11 and 10^'* erg 
for 11 < m < 100 in the second atep. Thia aaaiunes that 
10% of the eocpkiaioji energy of each SN (e.g., [6. S]} is used 
to drK'e the bulk, mot ton of the aiuToiiiiding gas. We talce 
{E^} = IC^ erg. Using the above results, we obtain 

,=s8.a.L48.)(^)(i^). im 

According to the theory of hierarchical atructuie forma^ 
tion in cold darft matter Dosmotogy, a Iialo with a total mass 
Mil coUapsing at redshift z > 1 has a radius 

whidi can be obtainal from Eq, 24 of [9] for the cmrent cos- 
moiogical parameters (fi^ = 0.27 and a TTubble constant of 

1. Kirfay EN. t^nliaiidi GA. SiiKin JD . Cnhcn JG. Guhsthalmrta P (2011) MtiU-ehHinnt 
abunilHKK En&acmnwrets from nKdun-nHokition 5p<tctr3. III. M«talEcTt]r dttrilHi- 
tions of MiUi) dwf saisUilc gatuiiH. Astinipti)! S 721-.TS. 

2. WiH J, CraiFtuu 5. DidEl A (2(108) Scdr^ idlalkin] jnd tho fumbnKnt^ he of ths 
loul dwf plaiiK. Men Mut Ro) AAicn SaL 1W:U53-1469. 

3. Mirth HF, dc Jong JTA, Ris M.W \2<Xa) A OHi^irshonrvi nuxbrun EkoElmil aful. 
JIOE of die itiuilinl praporticj of fasit MuBry Wjj saL^itAS. AxtfE4piTf i J 694:1075- 
1092. 

4. AndofT E. GrswiiH N (l^flS) AbuidanciH i^f tin alsmnti — Matati^k and soi^. 
Goodiin CiHiiwdiin: AcU S3:l 9T-214. 

£. W^OBfi S, NotiHlD K, Jaulu H-T, Kitam F&, Miilki B (IfflH] NuclHMjntllHHi Dn 
ohicinHi [3fitm njpaniovaG of asymfstotic .g^t bmdi ^tai. AjtfDp4Tji J 6l95:2Dd- 
220. 



TO km s ' Mpc '). Aa was necognEed in [10], thia purely 
theoretical result is in excellent agreement with what is re- 
quired (i.e., rh IX Afft'"^) to account for the slope of 3.S for the 
observed reSationship between log(M, /Af?.) and ([Fe/H]) for 
dSphs. We note that there is no a priori basis for dioosing an 
exphcil time (i.e.^ z) for the formation of dSphs. It is most 
plausible that they formed OEirlj in cold dark matter cosmol- 
ogy. Taking i + z = 10 and subatituting rh. from Eq, Sll into 
Eq. SIO. we obtain 

„ = 90.1Cl + 1.48t)(i^]'''. [12] 

In OUT mode!, {[Fe/H]) Iafi{Apb/[(l + fts 
logl0.aiei(l -I- lS.3fe)/Tj] (sec Eq. S7) and Af» fw 9.65 x 
W~'Mti/T} (see Eq. 14 in the main tact). As tj is deter- 
mined by Mh and k through Eq. SI 2. {[Ffe/H]) and M, are 
Functions of M^, and A'. We estimate and k Erom the data 
on ([Fe/H]) and Af, for eacb of the dSpha shown in Fig. 3 and 
give the nesulta in Tiible SI. Aa expected, these dSpha have 
very atrong SN-driven outflows with ij 10-360, The esti- 
mated values of AdT^ range from l.S x lO''^ A-fj, to 3.8 x 10* M^^ 
and are 3-70 times larger than the total mass within the 
haJf-lighL radius Afi/a inferred from observalionB [11]. The es- 
timated values of k range from 0.03 to 0.43 snd are typically 
-V 0.1. In particular, the relationship between log(A/, /AJq) 
and ({Ffe/H]) in Eq, S3 corresponds to *: = 0.14. 
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Table 1. Characteristics of dSphs 



Galaxies* 


([Fe/H])f 


(Afe) 


A^l/2* (Afg) 


Afh^ (A^a) 






inornax 


-0.99 


1.0 X 10'' 


T.4 X 10"" 


2.8- X 10" 


o.;30 


14 


Leo 1 


-1.43 


4.6 X 10* 


2.2 X lO'' 


1,1 X 10= 


0.14 


23 


Leo II 


-1.62 


1.4 X 10* 


7.3 X 10* 


5.4 X 10^ 


O.IB 


36 


Sculptor 


-1.68 


1.2 X 10* 


2.3 X 10'' 


4.S X 10^ 


0.13 


38 


Draco 


-1.93 


a.i X 10*^ 


2.1 X lO'' 


3.S X 10^ 


0.052 


40 


Sextans 


-1.93 


8.6 X 10* 


3.5 X lO'' 


3.6 X 10^ 


D.0E6 


41 


Una Minor 


-2.13 


s.e X 10* 


S.Sx lO'' 


2.S X 10^ 


0.023 


47 


Caies Venatici 1 


-1.98 


3.0 X 10* 


2.Sx 10'' 


2.0 X 10^ 


0.10 


SB 


Hercjles 


-2.41 


3.7 X 10* 


7.5 X 10* 


s.r X 10^ 


0.078 


147 


Ursa Major 1 


-2.18 


i.g X 10* 


1.3 X 10'' 


4,4 X 10^ 


0.31 


22B 


Leo IV 


-2.54 


8.5 X 10^ 


1.1 X 10* 


2.4 X 10^ 


0.13 


277 


Caies Venatici II 


-2.31 


8.0 X 10^ 


1.4 X 10* 


2,9 X 10^ 


D.4S 


351 


Urea Major II 


-2.47 


5.4 X 10^ 


7.9 X 10* 


2.0 X 10^ 


0.24 


357 


Coma Berenices 


-2.60 


4.8 X 10^ 


2.0 X 10* 


l.S X 10^ 


0.16 


354 



* Dmri sphatHdji j^iasiies in deoending oidercf M, . 
^Mean masllidiies ishsn Tram 

msmes uksn frDin |2] for ihe swen mosi m^Bsive dSphs and Tmni [3| fa ihe resi. 
^Tni^l dy[T.3nii[: rruBses wkhin h^lF-lighL radii uken from 
^EsiiniJKd uul masses of djH<: mjLTEr lialos. 
II Esiim jDed riljLivie rates Df SNe \i id CCENe. 
** Esiiiruied effidendes SN-driven auiAoNS. 
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